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This paper delineates the first steps in a systematic quantitative study of the spacetime fluctua- 
tions induced by quantum flelds in an evaporating black hole under the stochastic gravity program. 
The central object of interest is the noise kernel, which is the symmetrized two-point quantum cor- 
' relation function of the stress tensor operator. As a concrete example we apply it to the study of the 

, spherically-symmetric sector of metric perturbations around an evaporating black hole background 

• geometry. For macroscopic black holes we find that those fluctuations grow and eventually become 

04 ' important when considering sufficiently long periods of time (of the order of the evaporation time), 

but well before the Planckian regime is reached. In addition, the assumption of a simple correlation 
. between the fluctuations of the energy flux crossing the horizon and far from it, which was made 

' in earlier work on spherically-symmetric induced fluctuations, is carefully scrutinized and found to 

I be invalid. Our analysis suggests the existence of an infinite amplitude for the fluctuations when 

■ trying to localize the horizon as a three-dimensional hypersurface, as in the classical case, and, as 

' a consequence, a more accurate picture of the horizon as possessing a flnite effective width due 

to quantum fluctuations. This is supported by a systematic analysis of the noise kernel in curved 
spacetime smeared with different functions under different conditions, the details are collected in 
the appendices. This case study shows a pathway for probing quantum metric fluctuations near the 
horizon and understanding their physical meaning. 
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I. INTRODUCTION 



Studying the dynamics of quantum fields in a fixed curved spacetime, Hawking found that black holes emit thermal 
radiation with a temperature inversely proportional to their mass jlj . When the back reaction of the quantum fields 
I on the spacetime dynamics is included, one expects that the mass of the black hole decreases as thermal radiation 
■ at higher and higher temperatures is emitted. This picture, which constitutes the process known as black hole 
evaporation, is indeed obtained from semiclassical gravity calculations which are believed to be valid at least before 
the Planckian scale is reachedjl, Q . 



Semiclassical gravity [1, |^, [Jis a mean field description that neglects the fluctuations of the spacetime geometry. 
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, However, a number of studies have suggested the existence of large fluctuations near black hole horizons ^ ^ 

• (and even instabilities with characteristic time-scales much shorter than the black hole evaporation time. In all 
• ^ [ of them^ either states which are singular on the horizon (such as the Boulware vacuum for Schwarzschild spacetime) 
■ were explicitly considered, or fluctuations were computed with respect to those states and found to be large near the 
5h horizon. Whether these huge fluctuations are of a generic nature or an artifact arising from the consideration of states 
. . . I singular on the horizon is an issue that deserves further investigation. By contrast, the fluctuations for states regular 
on the horizon were estimated in Ref. [l^ and found to be small even when integrated over a time of the order of the 
evaporation time. 

These apparently contradictory claims and the fact that most claims on black hole horizon fluctuations were based 
on qualitative arguments and/or semi-quantitative estimates prompted us here to strive for a more quantitative and 
self-consistent description^. For this endeavor we follow the stochastic gravity program [l^ [l^ [13, [H, [ll] ■ We will 
consider the fluctuations of metric perturbations around a black hole geometry interacting with a quantum scalar field 
whose stress tensor drives the dynamics. The evolution of the mean background geometry is given by the semiclassical 
Einstein equation (with self-consistent back reaction from the expectation value of the stress tensor) while the metric 



^ At least those which provide a relativistic description. The argument in Refs. 0, H is based on a non-relativistic description and it 
is not obvious how to make some of our statements precise in that case. However, a natural generalization to the relativistic case is 
provided in Ref. [lOl . which does fall into this category. 

^ Previous attempts on this problem with similar emphasis by Raval, Sinha and one of us have appeared in Refs. Il3l |14|| . The apparent 
difference between the conclusions in Ref. [T3| and what is reported here will be explained below. 
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fluctuations obey an Einstein-Langevin equation [2l|, [S^, with a Gaussian stochastic source whose correlation 
function is given by the noise kernel, which characterizes the fluctuations of the stress tensor of the quantum fields. 
In contrast to the claims made before, we find here that even for states regular on the horizon the accumulated 
fluctuations become significant by the time the black hole mass has changed substantially, but well before reaching 
the Planckian regime. Our result is different from those obtained in prior studies, but in agreement with earlier work 
by Bekenstein ,24!]. 

The stochastic gravity program provides perhaps the best available framework to study quantum metric fluctuations, 
because while semiclassical gravity is a mean fleld description that does not take into account quantum metric 
fluctuations, the Einstein-Langevin equation enables one to solve for the dynamics of metric fluctuations induced by 
the fluctuations of the stress tensor of the quantum fields. Furthermore, the correlation functions that one obtains 
are equivalent to the quantum correlation functions for the metric perturbations around the semiclassical background 
that would follow from a quantum fleld theory treatment, up to a given order in an expansion in terms of the inverse 
number of fields [25l [26| . The quantization of these metric perturbations should be understood in the framework of 
a low-energy effective field theory approach to quantum gravity [l^l , which is expected to provide reliable results for 
phenomena involving typical length-scales much larger than the Planck length even if the microscopic details of the 
theory at Planckian scales are not known. 

A crucial relation assumed in previous investigations [l^ . [23 |'^ of the problem of metric fluctuations driven by 
quantum matter fleld fluctuations of states regular on the horizon (as far as the expectation value of the stress tensor 
is concerned) is the existence of correlations between the outgoing energy flux far from the horizon and a negative 
energy flux crossing the horizon, based on energy conservation arguments. Using semiclassical gravity, such correlations 
have been conflrmed for the expectation value of the energy fluxes, provided that the mass of the black hole is much 
larger than the Planck mass. However, a more careful analysis, summarized in Sec. IIVI shows that no such simple 
connection exists for energy flux fluctuations. It also reveals that the fluctuations on the horizon are in fact divergent 
unless it is treated as an object with a finite width rather than a three-dimensional hypersurface, as in the classical 
case, and one needs to find an appropriate way of probing the metric fiuctuations near the horizon and extracting 
physically meaningful information. This is a new challenge in the study of metric fluctuations which demands some 
careful thoughts on what they mean physically and how they can be probed operationally. In Appendices [B] and [C] 
we give a systematic analysis of the noise kernel in curved spacetime smeared with different functions under different 
conditions. The non-existence of this commonly invoked relation in this whole subject matter illustrates the limitations 
of heuristic arguments and the necessity of a detailed and consistent formalism to study the fluctuations near the 
horizon, in terms of their magnitude, how they are measured and their consequences. 

A few technical remarks are in order to delimit the problem under study: First, we will restrict our attention to the 
spherically-symmetric sector of metric fluctuations, which necessarily implies a partial description of the fluctuations. 
That is because, contrary to the case for semiclassical gravity solutions, even if one starts with spherically-symmetric 
initial conditions, the stress tensor fluctuations will induce fluctuations involving higher multipoles. Thus, the multi- 
pole structure of the fluctuations is far richer than that of spherically-symmetric semiclassical gravity solutions, but 
this also means that obtaining a complete solution (including all multipoles) for fluctuations rather than the mean 
value is much more difficult. 

Second, for black hole masses much larger than the Planck mass (otherwise the effective field theory description 
will break down anyway), one can introduce a useful adiabatic approximation involving inverse powers of the black 
hole mass. To obtain results to lowest order, it is sufficient to compute the expectation value of the stress tensor 
operator and its correlation functions in Schwarzschild spacetime. The corrections, proportional to higher powers of 
the inverse mass, can be neglected for sufficiently massive black holes. 

Third, when studying the dynamics of induced metric fiuctuations, the additional contribution to the stress tensor 
expectation value which results from evaluating it using the perturbed metric is often neglected. In the consideration 
of fiuctuations for an evaporating black hole such a term (which will be denoted by {Tj^l\g + h]) rcn in Sec. lIIip becomes 
important when it builds up for long times. The importance of this term is clear when comparing with the simple 
estimate made by Wu and Ford in Ref. [13, where {T^l\g + ^])rcn was neglected and the fluctuations were found to 
be small even when integrated over long times, of the order of the evaporation time of the black hole. 

The paper is organized as follows. In Sec.|TT]we briefly review the results for the evolution of the mean fleld geometry 
of an evaporating black hole obtained in the context of semiclassical gravity. The framework of stochastic gravity 



^ This approach has been mainly appHed to weak field situations, such as the study of quantum corrections to the Newtonian potential for 
particles in a Minkowski background [Ssl. l29l| . However, it is particularly interesting to apply it also to strong field situations involving 
cosmological [sOl or black hole spacetimes. 

See, however, Refs. [3ll. l32j. where those correlators were shown to vanish in an effectively two-dimensional model. 
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is then applied in Sec. IIIII to the study of the spherically-symmetric sector of fluctuations around the semiclassical 
gravity solution for an evaporating black hole. It has been previously assumed that an exact correlation between the 
fluctuations of the negative energy flux crossing the horizon and the flux far from it exists. In this paper we want to 
question this assumption, but in the presentation in Sec. IIIII we accept temporarily such a working hypothesis just 
so that we can have the common ground to compare our results with those in the literature. In Sec. IIVI we present 
a careful analysis of this assumption, and show that this supposition is invalid. Further details of this proof can be 
found in the Appendices [Bl and [Cl . Finally, in Sec.|V]we discuss several implications of our results and suggest some 
directions for further investigation. 

Throughout the paper we use Planckian units with h = c — G — 1 and the (+, +, +) convention of Ref. [3l|. We 
also make use of the abstract index notation of Ref. [13]. Latin indices denote abstract indices, whereas Greek indices 
are employed whenever a particular coordinate system is considered. 



II. MEAN EVOLUTION OF AN EVAPORATING BLACK HOLE 



Semiclassical gravity provides a mean field description of the dynamics of a classical spacetime where the gravita- 
tional back reaction of quantum matter fields is included self-consistently 0, H, @] . It is believed to be applicable to 
situations involving length-scales much larger than the Planck scale and for which the quantum back-reaction effects 
due to the metric itself can be neglected as compared to those due to the matter fields. The dynamics of the metric 
gab is governed by the semiclassical Einstein equation: 

Gab[9]^^(fab[g]) , (1) 
\ / ren 

where {Tab[g]) i-cn is the renormalized expectation value of the stress tensor operator of the quantum matter fields 
and K — STr/rrip with rrip being the Planck mass. One must solve both the semiclassical Einstein equation and the 
equation of motion for the matter fields evolving in that geometry, whose solution is needed to evaluate {Tab[g])ren 
self-consistently. 

An important application of semiclassical gravity is the study of black hole evaporation due to the back reaction 
of the Hawking radiation emitted by the black hole on the spacetime geometry. This has been studied in some detail 
for spherically symmetric black holes 0, ill . For a general spherically-symmetric metric there always exists a system 
of coordinates in which it takes the form 

ds^ = -e^^i^-^\l -2m{v, r) /r)dv^ + 2e'f'^'''''^ dvdr -f {dO^ + sin^ Odip^) . (2) 

This completely fixes the gauge freedom under local diffeomorphism transformations except for an arbitrary function 
of V that can be added to the function -(/)(«, r) and is related to the freedom in reparametrizing v (we will see below 
how this can also be fixed). In general this metric exhibits an apparent horizon, where the expansion of the outgoing 
radial null geodesies vanishes and which separates regions with positive and negative expansion for those geodesies, 
at those radii that correspond to (odd degree) zeroes of the vv metric component. We denote the location of the 
apparent horizon by rA_-a{v) — 2M{v), where M(v) satisfies the equation 2m{2M{v),v) = 2M{v). 

Spherical symmetry implies that the components Tgr, Tg^, T^r and T^t, vanish and the remaining components are 
independent of the angular coordinates. Under these conditions the various components of the semiclassical Einstein 
equation associated with the metric in Eq. ^ become 

Anr'T^, (3) 
-Wt;, (4) 
AirrTrr, (5) 

where in the above and henceforth we simply use T^j/ to denote the expectation value (r^iy[(7])ron and employ Planckian 
units (with — 1). Note that the arbitrariness in ip can be eliminated by choosing a parametrization of v such that 
tp takes a particular value at a given radius (we will choose that it vanishes at r = 2M(v), where the apparent horizon 
is located); is then entirely fixed by Eq. ([5]). 

Solving Eqs. ([SI)-® is no easy task. However, one can introduce a useful adiabatic approximation in the regime 
where the mass of the black hole is much larger than the Planck mass, which is in any case a necessary condition 
for the semiclassical treatment to be valid. What this entails is that when M ^ 1 (remember that we are using 
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Planckian units) for each value of v one can simply substitute T^i^ by its "parametric value" - by this we mean 
the expectation value of the stress energy tensor of the quantum field in a Schwarzschild black hole with a mass 
corresponding to M{v) evaluated at that value of v. This is in contrast to its dynamical value, which should be 
determined by solving self-consistently the semiclassical Einstein equation for the spacetime metric and the equations 
of motion for the quantum matter fields. This kind of approximation introduces errors of higher order in Ljj = B /M"^ 
{B is a dimensionless parameter that depends on the number of massless fields and their spins and accounts for their 
corresponding grey-body factors; it has been estimated to be of order 10~^ [SS^), which are very small for black holes 
well above Planckian scales. These errors are due to the fact that M{v) is not constant and that, even for a constant 
M(i)), the resulting static geometry is not exactly Schwarzschild because the vacuum polarization of the quantum 
fields gives rise to a non- vanishing {Tab[g\) ten [36]. 

The expectation value of the stress tensor for Schwarzschild spacetime has been found to cor resp ond to a thermal 
flux of radiation (with T:^ = Ln/i^Trr"^)) for large radii and of order Lr near the horizon^ [13, H US SO, EI ■ This 
shows the consistency of the adiabatic approximation for Lh ^ 1: the right-hand side of Eqs. I©-® contains terms 
of order Lh and higher, so that the derivatives of m(u, r) and ^{v, r) are indeed small. Furthermore, one can use the 
V component of the stress-energy conservation equation 

to relate the TJ" components on the horizon and far from it. Integrating Eq. ([6]) radially, one gets 

(r2r;)(r = 2M{v),v) - (r2T;)(r « 6Miv),v) + 0(L|), (7) 

where we considered a radius sufhciently far from the horizon, but not arbitrarily far (i.e. 2AI(v) <C r <C M{v)/Lii). 
The second condition is necessary to ensure that the size of the horizon has not changed much since the value of v' at 
which the radiation crossing the sphere of radius r at time v left the region close to the horizon. Note that while in the 
nearly flat region (for large radii) corresponds to minus the outgoing energy flux crossing the sphere of radius r, 
on the horizon, where ds^ = 2e'^^^'^''dvdr + {dO^ + s\v? 6d(p'^), equals Tyy, which corresponds to the null energy 
flux crossing the horizon. Hence, Eq. ([7]) relates the positive energy flux radiated away far from the horizon and the 
negative energy flux crossing the horizon. Taking into account this connection between energy fluxes and evaluating 
Eq. ^ on the apparent horizon, we flnally get the equation governing the evolution of its size: 

dM _ B 

Unless M(v) is constant, the event horizon and the apparent horizon do not coincide. However, in the adiabatic 
regime their radii are related, differing by a quantity of higher order in L^: r^niv) = rAuiv) (1 + ©(in))- 

We close this section with an explanation of why we did not have to deal with terms involving higher-order 
derivatives and even non-local terms when considering the expectation value of the stress tensor as one would expect 
for geometries with sufficiently arbitrary spacetime dependence of certain metric components. (This can be seen 
in explicit calculations for arbitrary Robertson- Walker geometries [lEI . l42j | or arbitrary small metric perturbations 
around specific backgrounds [il,|3|.) In our case such non-local and higher-order derivative terms would also appear 
in the exact expression of the stress tensor expectation value for arbitrary m(v,r) and ^(v,r) functions. However, 
the adiabatic approximation for M ^ 1 that we have employed effectively gets rid of them since one can replace 
the higher-derivative terms using Eqs. dSl)-® recursively and taking into account that the terms on the right-hand 
side are of order Therefore, higher-order derivative terms correspond to higher powers of 1/M^ and are 

highly suppressed for AI ^ 1. Note that this argument, which is based on the black hole size being much larger 
than the Planck length, is closely related to the order reduction prescription [0, that is often used to deal with 
higher-derivative terms in semiclassical gravity and other back-reaction problems. 



III. SPHERICALLY-SYMMETRIC INDUCED FLUCTUATIONS 



There are situations in which the fluctuations of the stress tensor operator and the metric fluctuations that they 
induce may be important, so that the mean fleld description provided by semiclassical gravity is incomplete and even 



^ The natural quantum state for a black hole formed by gravitational collapse is the Unruh vacuum, which corresponds to the absence 
of incoming radiation far from the horizon. The expectation value of the stress tensor operator for that state is finite on the future 
horizon of Schwarzschild, which is the relevant one when identifying a region of the Schwarzschild geometry with the spacetime outside 
the collapsing matter for a black hole formed by gravitational collapse. 
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fails to capture the most relevant phenomena (the generation of primordial cosmological perturbations constitutes a 
clear example of that). Stochastic gravity 1^ Ip . [itI Ip. lp| provides a framework to study those fluctuations. Its 



centerpiece is the Einstein-Langevin equation [20|, |21|, [22, 12. 

Gi'^ [g + h]^K (f^l^ [g + h]) +K u [g] , (9) 

\ / rcn 

which governs the dynamics of the metric fluctuations around a background metric gab that corresponds to a given 
solution of semiclassical gravity. The superindex (1) indicates that only the terms linear in the metric perturbations 
should be considered, and ^afc is a Gaussian stochastic source with vanishing expectation value and correlation function^ 
{iab{x)^cdix'))^ = {l/2){{iab{x),icdix')}) (with tab = Tab - (Tab)), whcrc the term on the right-hand side, which 
accounts for the stress tensor fluctuations within this Gaussian approximation, is commonly known as the noise 
kernel and denoted by Nabcd{x, x'). In this framework the metric perturbations are still classical but stochastic. 
Nevertheless, one can show that the correlation functions for the metric perturbations that one obtains in stochastic 
gravity are equivalent through order to the quantum correlation functions that would follow from a quantum 
field theory treatment when considering a large number of fields N p5l.[2^. In particular, the symmetrized two-point 
function consists of two contributions: intrinsic and induced fiuctuations. The intrinsic fiuctuations are a consequence 
of the quantum width of the initial state of the metric perturbations, and they are obtained in stochastic gravity by 
averaging over the initial conditions for the solutions of the homogeneous part of Eq. ([9]) distributed according to 
the reduced Wigner function associated with the initial quantum state of the metric perturbations. On the other 
hand, the induced fiuctuations are due to the quantum fluctuations of the matter fields interacting with the metric 
perturbations, and they are obtained by solving the Einstein-Langevin equation using a retarded propagator with 
vanishing initial conditions. 

In this section we study the spherically-symmetric sector [i.e., the monopole contribution, which corresponds to 
Z = 0, in a multipole expansion in terms of spherical harmonics y/„i(^, 0)] of metric fluctuations for an evaporating 
black hole. In this case only induced fluctuations are possible. The fact that intrinsic fluctuations cannot exist can be 
clearly seen if one neglects vacuum polarization effects, since Birkhoff 's theorem forbids the existence of spherically- 
symmetric free metric perturbations in the exterior vacuum region of a spherically-symmetric black hole that keep 
the ADM mass constant. Even when vacuum polarization effects are included, spherically-symmetric perturbations, 
characterized by m{v,r) and ip{v,r), are not independent degrees of freedom. This follows from Eqs. (l3])-([5]), which 
can be regarded as constraint equations. 

The fluctuations of the stress tensor are inhomogeneous and non-spherically-symmetric even if we choose a 
spherically-symmetric vacuum state for the matter fields (spherical symmetry simply implies that the angular de- 
pendence of the noise kernel in spherical coordinates is entirely given by the relative angle between the spacetime 
points X and x'). This means that, in contrast to the semiclassical gravity case, projecting onto the I = sector 
of metric perturbations does not give an exact solution of the Einstein-Langevin equation in the stochastic gravity 
approach that we have adopted here. Nevertheless, restricting to spherical symmetry in this way gives more accurate 
results than two-dimensional dilaton-gravity models resulting from simple dimensional reduction [46l . |47| . |48| . This 
is because we project the solutions of the Einstein-Langevin equation just at the end, rather than considering only 
the contribution of the s-wave modes to the classical action for both the metric and the matter fields from the very 
beginning. Hence, an infinite number of modes for the matter fields with / ^ contribute to the I — projection of the 
noise kernel, whereas only the s-wave modes for each matter field would contribute to the noise kernel if dimensional 
reduction had been imposed right from the start, as done in Refs. [3l|,[32,|4^ as well as in studies of two-dimensional 
dilaton-gravity models. 

The Einstein-Langevin equation for the spherically-symmetric sector of metric perturbations can be obtained by 
considering linear perturbations of m{v,r) and ip{v,r), projecting the stochastic source that accounts for the stress 
tensor fiuctuations to the I = sector, and adding it to the right-hand side of Eqs. ©-([S]). We will focus our attention 
on the equation for the evolution of 77(1;, r), the perturbation of m{v,r): 

which reduces, after neglecting terms of order or higher, to the following equation to linear order in r]: 

|!^ = ^, + Wc (11) 



Throughout the paper we use the notation (...){ for stochastic averages over all possible realizations of the source ^^i) to distinguish 
them from quantum averages, which are denoted by (...). 
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It is important to emphasize that in Eq. (jTO]) we assumed that the change in time of ri{v, r) is sufficiently slow so 
that the adiabatic approximation employed in the previous section to obtain the mean evolution of m{v,r) can also 
be appHed to the perturbed quantity m{v,r) + r]{v,r). This is guaranteed as long as the term corresponding to the 
stochastic source is of order Lh or higher, a point that will be discussed below. 

Obtaining the noise kernel which determines the correlation function for the stochastic source is highly nontrivial 
even if we compute it on the Schwarzschild spacetime, which is just ified in the adiabatic regime for the background 
geometry. As implicitly done in prior work (for instance in Refs. |12|.[2^: see, however, Refs. [Sllls^). we will assume 
in this section that the fluctuations of the radiated energy flux far from the horizon are exactly correlated with the 
fluctuations of the negative energy flux crossing the horizon. This is a crucial assumption which implies an enormous 
simplification and allows a direct comparison with the results in the existing literature, and its validity will be analyzed 
more carefully in the next section.^ 

Since the generation of Hawking radiation is especially sensitive to what happens near the horizon, from now on 
we will concentrate on the metric perturbations near the horizon^ and consider 'r]{v) = r](v,2M(v)). Assuming that 
the fluctuations of the energy flux crossing the horizon and those far from it are exactly correlated, from Eq. (fTTj) we 
have 

where ^(w) = (47rr^ r « 6M{v)). The correlation function for the spherically-symmetric fluctuation is 

determined by the integral over the whole solid angle of the iV'^'^ component of the noise kernel, which is given by 
{l/2){{t'^{x),tl{x')}). The I = fluctuations of the energy flux of Hawking radiation far from a black hole formed 
by gravitational collapse, characterized also by {l/2){{t'^{x),t^{x')}) averaged over the whole solid angle, have been 
studied in Ref. [12]. Its main features are a correlation time of order M and a characteristic fluctuation amplitude 
of order eo/M^ (this is the result of smearing the stress tensor two-point function, which diverges in the coincidence 
limit, over a period of time of the order of the correlation time). The order of magnitude of eg has been estimated to 
lie between O.IB and B [13, [HI- For simplicity, we will consider quantities smeared over a time of order M. We can 
then introduce the Markovian approximation {eo/M^{v))S{v — v'), which coarse-grains the information on features 
corresponding to time-scales shorter than the correlation time M . Under those conditions r^^^' is of order and 
thus the adiabatic approximation made when deriving Eq. (I10|) is justified. 

The stochastic equation p2p can be solved in the usual way and the correlation function for ri{v) can then be 
computed. Alternatively, one can follow Bekenstein [23 | and derive directly an equation for (77^ (u))^. This is easily 
done multiplying Eq. p2p by ry(u) and taking the expectation value. The result is 

|;(^' W>« = + mvnv))^. (13) 

For delta-correlated noise (the Stratonovich prescription is the appropriate one in this case), (?7('u)^(w))^ equals one 
half the time-dependent coefhcient multiplying the delta function 5(v — v') in the expression for {S.iv)^{v'))^, which 
is given by eo/M^{v) in our case. Taking that into account, Eq. becomes 

Finally, it is convenient to change from the v coordinate to the mass function M{v) for the background solution. 
Eq. (HH) can then be rewritten as 

_^(,„.,M,,,.-A(,»,M„,-Ma^ (15. 



This simple relation between the energy flux crossing the horizon and the flux far from it is valid for the expectation value of the stress 
tensor, which is based on an energy conservation argument for the component. In most of the literature this relation is assumed to 
hold also for fluctuations. However, in the next section we will show that this is an incorrect assumption. Therefore, results derived 
from this assumption and conclusions drawn are in principle suspect. (This misstep is understandable because most authors have not 
acquired as much insight into the nature of fluctuations phenomena as now.) Our investigation testifies to the necessity of a complete 
reexamination of all cases afresh. In fact, an evaluation of the noise kernel near the horizon seems unavoidable for the consideration of 
fluctuations and back-reaction issues. 
* This means that possible effects on the Hawking radiation due to the fluctuations of the potential barrier for the radial mode functions 
will be missed by our analysis. 
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The solutions of this equation are given by 



M AB 



Mo 
M 



(16) 



Provided that the fluctuations at the initial time corresponding to M — Mq are negligible (much smaller than 
^/ Co /4:B ~ 1), the fluctuations become comparable to the background solution when M ~ Mq^^. Note that fluctu- 
ations of the horizon radius of order one in Planckian units do not correspond to Planck scale physics because near 
the horizon AR — r — 2M corresponds to a physical distance L ^ V M AR, as can be seen from the line element 
for Schwarzschild, ds^ — —(1 — 2M/r)dt^ + (1 ^ 2M/r)^^dr'^ + r'^{d9'^ + sir? 6 d^p^), by considering pairs of points 
at constant t. So AR ~ 1 corresponds to L ~ whereas a physical distance of order one is associated with 

AR ~ 1/M, which corresponds to an area change of order one for spheres with those radii. One can, therefore, have 
initial fluctuations of the horizon radius of order one for physical distances well above the Planck length provided that 
we consider a black hole with a mass much larger than the Planck mass. One expects that the fluctuations for states 
that are regular on the horizon correspond to physical distances not much larger than the Planck length, so that the 
horizon radius fluctuations would be much smaller than one for sufflciently large black hole masses. Nevertheless, that 
may not be the case when dealing with states which are singular on the horizon, with estimated fluctuations of order 
Qj. even y/M [HI- Confirming that the fluctuations are indeed so small for regular states and verifying 

how generic, natural and stable they are as compared to singular ones is a topic that we plan to address in future 
investigations. 

Our result for the growth of the fluctuations of the size of the black hole horizon agrees with the result obtained 
by Bekenstein in Ref. [2^ and implies that, for a sufficiently massive black hole (with a few solar masses or a 
supermassive black hole), the fluctuations become important before the Planckian regime is reached. Strictly speaking, 
one cannot expect that a linear treatment of the perturbations provides an accurate result when the fluctuations 
become comparable to the mean value, but it signals a significant growth of the fluctuations (at least until the 
nonlinear effects on the perturbation dynamics become relevant). 

This growth of the fluctuations which was found by Bekenstein and conflrmed here via the Einstein-Langevin 
equation seems to be in conflict with the estimate given by Wu and Ford in Ref. [13] ■ According to their estimate, 
the accumulated mass fluctuations over a period of the order of the black hole evaporation time {At ~ Mg ) would 
be of the order of the Planck mass. The discrepancy is due to the fact that the first term on the right-hand side of 
Eq. which corresponds to the perturbed expectation value {t'^ [g + h]) in Eq. was not taken into account 
in Ref. [I4I . The larger growth obtained here is a consequence of the secular effect of that term, which builds up in 
time (slowly at first, during most of the evaporation time, and becoming more significant at late times when the mass 
has changed substantially) and reflects the unstable nature of the background solution for an evaporating black hole.^ 

All this can be qualitatively understood as follows. Consider an evaporating black hole with initial mass Mq and 
suppose that the initial mass is perturbed by an amount 5Mq = 1. The mean evolution for the perturbed black hole 
(without taking into account any fluctuations) leads to a mass perturbation that grows like 5M = {Mq/M)'^ SMq = 
{Mq/M)"^ , so that it becomes comparable to the unperturbed mass M when M ~ A/q^^, which coincides with the result 
obtained above. Such a coincidence has a simple explanation: the fluctuations of the Hawking flux slowly accumulated 
during most of the evaporating time, which are of the order of the Planck mass, as found by Wu and Ford, give a 
dispersion of that order for the mass distribution at the time when the instability of the small perturbations around 
the background solution start to become signiflcant. 



IV. CORRELATION BETWEEN OUTGOING AND INGOING ENERGY FLUXES 



In this section we address the issue whether the simple relation between the energy flux crossing the horizon and 
the flux far from it also holds for the fluctuations. As we show in Appendix[Xl those correlations vanish for conformal 
fields in two-dimensional spacetimes. (The correlation function for the outgoing and ingoing null energy fluxes in an 
effectively two-dimensional model was explicitly computed in Refs. [3ll,[32| and it was indeed found to vanish.) On the 



^ A clarification between our results and the claims by Hu, Raval and Sinha in Ref. is in place here: both use the stochastic gravity 
framework and perform an analysis based on the Einstein-Langevin equation, so there should be no discrepancy. However, the claim 
in Ref. [l3l | was based on a qualitative argument that focused on the stochastic source, but missed the fact that the perturbations 
around the mean are unstable for an evaporating black hole. Once this is taken into account, agreement with the result obtained here 
is recovered. 
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other hand, in four dimensions the correlation function does not vanish in general and correlations between outgoing 
and ingoing fluxes do exist near the horizon (at least partially). This point is also explained in Appendix \X\ 

For black hole masses much larger than the Planck mass, one can use the adiabatic approximation for the background 
mean evolution. Therefore, to lowest order in Ljj one can compute the fluctuations of the stress tensor in Schwarzschild 
spacetime. In Schwarzschild, the amplitude of the fluctuations of r^TJ" far from the horizon is of order 1/AP {— 
JVP/M'^) when smearing over a correlation time of order M, which one can estimate for a hot thermal plasma in flat 
space [13, IMi (see also Ref. for a computation of the fluctuations of r^T^ far from the horizon). The amplitude of 
the fluctuations of r^T^ is thus of the same order as its expectation value. However, their derivatives with respect to 
V are rather different: since the characteristic variation times for the expectation value and the fluctuations are 
and M respectively, d{r^T^)/dv is of order whereas 9(r^^JJ)/9u is of order This implies an additional 

contribution of order L^i due to the second term in Eq. ^ if one radially integrates the same equation applied to 
stress tensor fluctuations (the stochastic source in the Einstein-Langevin equation). Hence, in contrast to the case 
of the mean value, the contribution from the second term in Eq. ([5]) cannot be neglected when radially integrating 
since it is of the same order as the contributions from the flrst term, and one can no longer obtain a simple relation 
between the outgoing energy flux far from the horizon and the energy flux crossing the horizon. 

So far we have argued that the method employed for the mean value cannot be employed for the fluctuations. 
Although one expects that r^^JJ on the horizon and far from it will not be equal when including the contributions 
that results from radially integrating the second term in Eq. ([6]), one might wonder whether there is a possibility 
that those contributions would somehow cancel out. That possibility can, however, be excluded using the following 
argument. The smeared correlation function 

J dvh{v) J dv'h{v')r^{C{v,r)C{v',r))^, (17) 

where h{v) is some appropriate smearing function and ^y{v, r) has already been integrated over the whole solid angle, 
is divergent on the horizon but flnite far from it. Therefore, r^^^ on the horizon and far from it cannot be equal for 
each value of v. 

Let us discuss in some more detail the fact that certain smearings of the quantity r'^{£,l{v, r)^l{v' , r))j are divergent 
on the horizon but flnite far from it. The smeared correlation function is related to the noise kernel as follows: 

j dvdv'h{v)h{v') r^{CAv, r)ilW. r))(_ = J dvdv' h{v)h{v') J dndfl'NY^iv, r, 6, ^; v' , r, 0', ip'). (18) 

The noise kernel is divergent in the coincident limit or for null-separated points. Smearing the noise kernel along all 
directions gives a flnite result. However, although certain partial smearings also give a flnite result, others do not. For 
instance, smearing along a timelike direction yields a finite result, whereas smearing on a spacelike hypersurface yields 
in general a divergent result [s^]- On the other hand, the result of smearing along two "transverse" null directions 
(two null directions sharing the same orthogonal spacelike 2-surfaces) is also flnite, but not for a smearing along just 
one null direction even if we also smear along the orthogonal spacelike directions. For r > 2M Eq. corresponds to 
a smearing along a timelike direction and gives a finite result for the smeared correlation function, but on the horizon 
it corresponds to a smearing along a single null direction and it is divergent. 

A proof of the results described in the previous paragraph is provided in Appendices [B] and [C] by considering a 
product of smearing functions involving all directions and then taking different kinds of limits in which the smearing 
size along certain directions vanishes. One limit corresponds to taking a vanishing size for the smearing along one of 
the null directions and we refer to it as smearing along null geodesies. The other corresponds to taking vanishing sizes 
for certain spatial directions and we refer to it as smearing along timelike curves. The proof proceeds in two steps. 
First, in Appendix [Bl it is shown for the fiat space case. Then it is generalized to curved spacetimes in Appendix [Cl 
using a quasilocal expansion in terms of Riemann normal coordinates. 

V. DISCUSSION 

Following the stochastic gravity program, in Sec. IHII we found that the spherically-symmetric fiuctuations of the 
horizon size of an evaporating black hole become important at late times, and even comparable to its mean value when 

2/3 

M ^ Mq , where Mq is the mass of the black hole at some initial time when the fluctuations of the horizon radius 
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are much smaller than the Planck length. -"^^ This is consistent with the result previously obtained by Bekenstein in 



It is important to realize that for a sufficiently massive black hole, the fluctuations become significant well before 
the Planckian regime is reached. More specifically, for a solar mass black hole they become comparable to the mean 
value when the black hole radius is of the order of lOnm, whereas for a supermassive black hole with M ^ 10^ Mq, that 
happens when the radius reaches a size of the order of 1mm. One expects that in those circumstances the low-energy 
effective field theory approach of stochastic gravity should provide a reliable description. 

It is worth mentioning that other properties of the black hole can exhibit substantial fluctuations when taking into 
account the back reaction of Hawking radiation. As pointed out by Page in Ref. [ssj, by momentum conservation the 
fluctuations of the total momentum of the Hawking radiation emitted will cause a recoil of the position of the black 
hole, which will also fluctuate. According to Page's estimate, the spread of the distribution for the black hole position 
will become comparable to the size of the horizon by the time an energy of order M^/'^ has been emitted. This kind of 
fluctuations, which were not obtained in our calculation because we restricted our attention to spherically-symmetric 
metric perturbations, exhibit certain features that differ significantly from those of our result. The fluctuations that 
we found take a much longer time to build up and depend crucially on the unstable behavior of small perturbations of 
the semiclassical solution, with characteristic time-scales of the order of the evaporation time. On the contrary, this 
unstable behavior plays no major role in the growth of the position fluctuations. Furthermore, since this growth is still 
much slower than the emission rate of Hawking quanta, in the frame where the black hole is at rest the properties of the 
Hawking radiation being emitted remain essentially unchanged when the position fluctuations are taken into account, 
whereas the fluctuations of the horizon size do imply fluctuations in the temperature of the Hawking radiation. 

Due to the nonlinear nature of the back-reaction equations, such as Eq. pi}|) . the fact that the fluctuations of the 
horizon size can grow and become comparable to the mean value implies non-negligible corrections to the dynamics 
of the mean value itself. This can be seen by expanding Eq. (jlOp (evaluated on the horizon) in powers of rj and taking 
the expectation value. Through order rf we get 



When the fluctuations become comparable to the mass itself, the third term (and higher order terms) on the right- 
hand side is no longer negligible and we get non-trivial corrections to Eq. ([5]) for the dynamics of the mean value. 
These corrections can be interpreted as higher order radiative corrections to semiclassical gravity that include the 
effects of metric fluctuations on the evolution of the mean value. For instance, the third term on the right-hand side 
of Eq. (jl9p would correspond to a two-loop Feynman diagram involving a matter loop with an internal propagator for 
the metric perturbations (restricted to the spherically-symmetric sector in our case), as compared to just one matter 
loop, which is all that semiclassical gravity can account for. 

Does the existence of the significant deviations for the mean evolution mentioned above imply that the results 
based on semiclassical gravity obtained by Bardeen and Massar in Refs. [2j l2J are invalid? Several remarks are in 
order. First of all, those deviations start to become significant only after a period of the order of the evaporation 
time when the mass of the black hole has decreased substantially. Secondly, since fluctuations were not considered in 
those references, a direct comparison cannot be established. However, we can compare the average of the fluctuating 
ensemble with their results. Doing so exhibits an evolution that deviates signiflcantly when the fluctuations become 
important. Nevertheless, if one considers a single member of the ensemble at that time, its evolution will be accurately 
described by the corresponding semiclassical gravity solution until the fluctuations around that particular solution 
become important again, after a period of the order of the evaporation time associated with the new initial value of 
the mass at that time. 

An interesting aspect that we have not addressed in this work, but which is worth investigating, is the quantum 
coherence of those fluctuations. It seems likely that, given the long time periods involved and the size of the fluctua- 
tions, the entanglement between the Hawking radiation emitted and the black hole spacetime geometry will effectively 
decohere the large horizon fluctuations, rendering them equivalent to an incoherent statistical ensemble. 

In this paper we have taken a first step to put the study of metric fluctuations in black hole spacetimes on a firmer 
basis by considering a detailed derivation of the results from an appropriate formalism rather than using heuristic 



Remember that for large black hole masses this can still correspond to physical distances much larger than the Planck length, as 
explained in Sec. IIIII 



Ref. [2^. 




(19) 
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arguments or simple estimates. The spirit is somewhat analogous to the study of the mean back-reaction effect of 
Hawking radiation on a black hole spacetime geometry (both for black holes in equilibrium and for evaporating ones) 
by considering the solutions of semiclassical gravity in that case rather than just relying on simple energy conservation 
arguments. In order to obtain an explicit result from the stochastic gravity approach and compare with earlier work, in 
Sec, mil we employed a simplifying assumption implicitly made in most prior work: the existence of a simple connection 
between the outgoing energy flux fluctuations far from the horizon and the negative energy flux fluctuations crossing 
the horizon. In Sec. IIVI we analyzed this assumption carefully and showed it to be invalid. This strongly suggests 
that one needs to study the stress tensor fluctuations from an explicit calculation of the noise kernel near the horizon. 
This quantity is obtainable from the stochastic gravity program and calculation is underway jUj, [H, [5^ . 

A possible way to compute the noise kernel near the horizon could be to use an approximation scheme based on 
a quasilocal expansion such as Page's approximation [38] or similar methods corresponding to higher order WKB 
expansions [4l|.^^ With these techniques one can obtain an approximate expression for the Wightman function of the 
matter fields, which is the essential object needed to compute the noise kernel. Unfortunately these approximations 
are only accurate for pairs of points with a small separation scale and break down when it becomes comparable to 
the black hole radius. Therefore, especial care is needed when studying the I = multipole since that corresponds 
to averaging the noise kernel over the whole solid angle, which involves typical separations for pairs of points on 
the horizon of the order of the black hole radius, and one needs to make sure that the integral is dominated by the 
contribution from small angular separations. Alternatively, one might hope to gain some insight on the fluctuations 
near a black hole horizon by studying the fluctuations of the event horizon surrounding any geodesic observer in de 
Sitter spacetime, which exhibits a number of similarities with the event horizon of a black hole in equilibrium [STj . 
In contrast to the black hole case, it may be possible to obtain exact analytical results for de Sitter space due to its 
high degree of symmetry. 

Furthermore, as explained in Sec. IIVI and shown in detail in Appendices |B] and [Cl the noise kernel smeared over 
the horizon is divergent, and so are the induced metric fluctuations. Hence, one cannot study the fluctuations of 
the horizon as a three-dimensional hypersurface for each realization of the stochastic source because the amplitude 
of the fluctuations is infinite, even when restricting one's attention to the I = sector. Instead one should regard 
the horizon as possessing a finite effective width due to quantum fluctuations. In order to characterize its width one 
must find a sensible way of probing the metric fluctuations near the horizon and extracting physically meaningful 
information, such as their effect on the Hawking radiation emitted by the black hole. One possibility is to study how 
metric fluctuations affect the propagation of a bundle of null geodesies [31, 32, 49, [s^ [H, '60]. One expects that 
this should provide a description of the effects on the propagation of a test field whenever the geometrical optics 
approximation is valid. However, preliminary analysis of simpler cases with a quantum field theory treatment suggest 
that when including quantum vacuum metric fluctuations the geometrical optics approximation becomes invalid. 
Another possibility, which seems to constitute a better probe of the metric fluctuations, is to analyze their effect on 
the two-point quantum correlation functions of a test field. The two-point functions characterize the response of a 
particle detector for that field and can be used to obtain the expectation value and the fluctuations of the stress tensor 
of the test fleld. 

Finally, since the large fluctuations suggested in Refs. 0, H, S [13] involve time-scales much shorter than the 
evaporation time (contrary to those considered in this paper) and high multipoles, one expects that for a sufficiently 
massive black hole the spacetime near the horizon can be approximated by Rindler spacetime (identifying the black 
hole horizon and the Rindler horizon) provided that we restrict ourselves to sufficiently small angular scales. Thus, 
analyzing the effect of including the interaction with the metric fluctuations on the two-point functions of a test field 
propagating in flat space, which is technically much simpler, could provide useful information for the black hole case. 
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Note, however, that in most of these approaches the state of the quantum fields is the Hartle-Hawking vacuum. For an evaporating 
black hole, one should consider the Unruh vacuum. 
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APPENDIX A: CORRELATIONS BETWEEN OUTGOING AND INGOING FLUXES IN 1 + 1 AND 3 + 1 

DIMENSIONS 

Any spacetime metric in 1 + 1 dimensions is conformally flat (at least locafly) and can be written as 

ds^ ^ ~C{u,v)dudv. (Al) 
In terms of these null coordinates, the conservation equation for the stress tensor, V°'Tab — 0, reduces to 

since all the other relevant Christoffel symbols vanish. Taking into account that r„„ — 9„(lnC), F^^ — dv{\nC) and 
combining Eqs. (IX2|l - (IX3|) we get 

a„f„„ + 9„f,„ = ~Cdt{f^.JC), (A4) 
dvTuu ~ duT^v = C dx {Tuv / C) , (^5) 

where we introduced the coordinates t = {u + v)/2 and x = {v — u)/2. 

This result can be applied to the Schwarzschild geometry in 1 + 1 dimensions identifying t with the usual Killing 
time and x with the Regge- Wheeler coordinate r*. In that case we have C{t,r) = (1 — 2M/r). If we consider a 
massless conformally coupled field in 1 + 1 dimensions (conformal and minimal coupling are equivalent in that case) , 
the trace of the stress tensor (which is related to the Tuv component) vanishes at the classical level and is entirely 
given by the trace anomaly {T^)i-cn — R/24tt = M/Girr^. Since both the trace anomaly and the conformal factor 
are time-independent, the term on the right-hand side of Eqs. (jA4p vanishes. Therefore, it follows from Eq. (|A4p 
that the generation of left-moving and right-moving null mean fluxes is perfectly anticorrelated, which implies that 
the positive energy flux of outgoing Hawking radiation equals, in absolute value, the negative energy flux crossing 
the horizon. Moreover, from Eq. (|A4|) and the fact that {Tuv)rcn/C = — (r^)rcn/4 — — Af/247rr'^, which implies 

C drX{Tuv)i-cn/C) = M/Qnr^, it is clear that the amount of anticorrelated mean fluxes generated tends to zero for 
large radii and is largest at r = 3M. 

For energy flux fluctuations, the situation is very different. The trace anomaly does not fluctuate, i.e., the trace 
of the stress tensor does not fluctuate for conformal flelds [l^. Hence, in 1 -I- 1 dimensions neither correlated nor 
anticorrelated fluctuations of the left-moving and right-moving null fluxes can be generated in the absence of other 
interactions: both fluxes are separately conserved. 

On the other hand, although in 3-1-1 dimensions one can try to use a similar argument when considering sectors with 
certain symmetries, the final conclusion is different. For instance, if one uses in Minkowski spacetime a coordinate 
system {u,v,y,z\ where u = t — x and v — t + x are null coordinates, the components Tuy, T^z, T^y, Ty^ and Tyz 
vanish in the sector which is rotationally invariant on the yz plane. In that sector, the u and v components of the 
conservation equation coincide with those in the 1-1-1 case, given by Eqs. (|A2[) - (jA3p . with vanishing Christoffel symbols. 
Therefore, the conditions for the generation of correlated or anticorrelated null fluxes are still given by Eqs. (jA4[) - (jA5p . 
However, in 3 -I- 1 dimensions the Tyy and Tzz components also contribute to the trace of the stress tensor. Hence, 
although the trace does not fluctuate for conformal fields, Tuv does fluctuate and its fluctuations coincide with those 
of {Tyy + Tzz)/'^- Correlated and anticorrelated null energy flux fluctuations can thus be generated. 

The discussion in the previous paragraph can be extended to a general spherically symmetric spacetime in the 
region where the gradient of the radial coordinate is spacelike. This can be done by considering angular coordinates 
for every sphere of constant radius and the coordinates associated with the two radial null directions orthogonal to 
them. Taking rotational symmetry into account, as done in the flat space case, one can provide an argument similar 
to that in 1 + 1 dimensions. However, as in the flat space case, despite the absence of fluctuations in the trace of the 
stress tensor, T^v will still fluctuate due to the fluctuations of {Tgg+T^^/ sin^ 9)/'ir^. Moreover, in this case there will 
be additional contributions due to scattering off the potential barrier (an effect that will be small near the horizon). 

APPENDIX B: SMEARING OF THE NOISE KERNEL IN FLAT SPACE 

In this appendix we will consider several kinds of smearings of the noise kernel for the Minkowski vacuum. In Sec. IB II 
we study a product of smearing functions involving two null directions and the two orthogonal spatial directions, and 
analyze the limit in which the smearing size along one of the null directions vanishes, which is shown to be divergent. 
On the other hand, in Sec. IB 21 a smearing along a timelike direction and the three orthogonal directions is considered 



F" T — n 

^ vv vv 



(A2) 
(A3) 
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and it is shown that a finite smearing size along the timehke direction is sufficient to have a finite result. In contrast, in 
the limit of a vanishing smearing size along the timelike direction the result is always divergent, even for non- vanishing 
smearing sizes along all the spatial directions. 

The noise kernel should be treated as distribution in spacetime coordinates. It has a divergent coincidence limit 
and it involves subtle integration prescriptions. However, the Fourier transforms of this kind of distributions are much 
simpler to deal with. Therefore, it is very convenient to perform our calculations in Fourier space and we will do so 
throughout this appendix. 

The noise kernel Nabcd{x,x') — {l/2){{Tab{x),Tcd{x')}) — {Tab{x)){Tcd{x')) for a massless conformally coupled 
scalar field in the Minkowski vacuum state has been obtained for instance in Ref. [61] and is given in standard inertial 
coordinates by 

N^,pa{x, x') = ^ (3I?^(pI?,), - V^.Vp,) N{x - x'), (Bl) 

where V^^, = -q^yU - d^^d^ and 

A^(-x')=(^/|fje-(-').(V). (B2) 



1. Smearing around a null geodesic 

In this subsection we consider the case in which one smears around a null geodesic. Using the null coordinates 
V — t + X and u = t — x, we define the smeared version of the kernel N{x — x') as 

7V=y dudvdPxg{u)h{v)f{x) j du'dv'd^x'g{u')h{v')f{x')N{v-v',u-u',x~x'), (B3) 

where we integrated the kernel with some smearing functions both for the null coordinates u and v, and for the 
orthogonal spatial directions. If we choose Gaussian smearing functions 

giu) = {2TT<jlr-2 exp{-uy2al), (B4) 

h{v) = {2Tral)--2 eM-vV2al), (B5) 

fix) = (2^a,2)-iexp(-fV2a2), (B6) 



Eq. (IB3P can be written as 

= ,4 \ 2 4 j dYdUd^Xe'^e'^ke'^ / dA.dA.^d^A e"^e"^e"^7V(A„ A„, A), (B7) 
(27r)4cr^;cr,>4 J J 

where we introduced the semisum and difference variables U ~ {u + u')/2, V — {v + v')/2, X ~ {x + x')/2, Au — u' — u, 
Ay = v' ~ V and A = x' — x. The integrals for [/, V and X can be readily performed and yield the following result: 

1 r _ _ _ 

^=7T^2 ^ dAydAud^Ae ^iV(A,, A„, A). (B8) 

On the other hand, in order to compute the remaining integrals it is convenient to work in Fourier space. When 
considering the null coordinates v and u, it is useful to introduce the momenta py — {px —pt)/2 and p„ = —{pt +Px)/2 
so that Eq. (|B2|1 becomes 



where we used the vector notation for the transverse components associated with the coordinates y and z. Eq. (jB8[) 
can then be expressed in Fourier space as 

A/" = (1920^) (2^)4 / dp«rfp«rf'pe-''^'^"e"P"-"e-P"'^?- - f). (BIO) 
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One can then infer that M diverges as (Ju 0. This can be seen as follows. The integral in Eq. (jBlOp gets two 
identical contributions from the quadrants ipu,Pv > 0) and (pu,Pv < 0), whereas the remaining two quadrants give 
a vanishing contribution. Moreover, since the integrand is positive, the integral is also positive and greater than the 
same integral restricted to a smaller domain of integration. Taking all that into account, we have 



PvPu >a-r 
Pv - Pu >0 



2e^i 1 f°° 2 2 Z"'^" 

(19207r)(27r)4 72^^,/^2 J,y-'/2 
1 r°° 2 2 1 

= 2 noon V9 H / -. (Bll) 

cr^CTi, (19207r)(27r)4 _/2<T„/<x2 (Juaya^ 

The last integral is divergent if one takes (t„ ^ (at least for <Tr 7^ ^'^). Thus, Af diverges unless (7„ 7^ 0. 

Using the previous result, it is easy to discuss whether a smeared version of the actual noise kernel (including the 
differential operators), given by Eq. (|Bip . also diverges when i7„ — > 0. Each derivative in Eq. (jBl[) gives rise to an 
additional factor involving the momentum associated with the corresponding component. Additional factors involving 
powers of p„, p^, Py and p^ (odd powers of Py or p^ give a vanishing contribution) leave the argument employed in 
the previous paragraph unchanged and the same conclusions obtained when a„ — > hold for the smeared noise kernel 
as well (including the differential operators). For example we have Afyvw i/<^uO'yCr'^ for the smeared version of the 
WW component of the noise kernel. 

On the other hand, when both ay and (Tu 7^ the smeared kernel Af is finite even for ar = 0. This can be seen 
by taking (7^ = in Eq. (jBlOp . We then have 



^ = (1920.) (2^)4 / dpM'pe-P'^'e-^'^ieipyPy-p') 
(ipt,dp„7rp„p„e~P"'^-e~^'"'^" 0{PvPu) 



1 1 



(19207r) (27r)4 
2 1 



(1920.) (27r)4 7o 



dpy I dp„7rp„p„e-''>"e-f"'^' ~ (B12) 



which is finite for au,crv 7^ 0. It is also clear that the same conclusion applies to the smeared version of the noise 
kernel. For instance we have Mww ^ l/^u'^^j- 

2. Smearing along a timelike curve and on a spacelike hypersurface 

In this subsection we consider the smeared kernel M obtained when working with standard cartesian coordinates 
in Minkowski spacetime: 

M= I dt(fxg{t)f{x) [ dt'd^x'g{t')f{x')N{t-t\x- x'), (B13) 



where we used the vector notation for the three spatial components a;, y and z, and the Gaussian smearing functions 
are now 

g{t) = (27ra2)-i exp(-iV2a?), (B14) 
fix) = (27ra2)-i exp(-fV2^T')- (B15) 

Note that Ur corresponds now to the smearing size of the three spatial directions. Introducing the semisum and 
difference variables T = [t + t')/2 and At =t' — t, Eq. (|B13[) can be rewritten as 

1 r f _ a2 

A/"= ,^ ,, „ ^ dTd^Xe ^ dAtd^Ae ^iV(At,A). (B16) 

(27r)''crj o-o J J 



The divergence of as (t^ — > can also be proven for <Jr = {i.e., in the absence of smearing along the transverse spatial directions) 
by taking Cj. = in Eq. IIBlOl l and replacing crj; with an arbitrary but fixed positive value in Eq. IIBlll l. 
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After integrating over T and X we have 



■^^TT^I T I dAtd^Ae ^7V(At,A), (B17) 



which can be equivalently rewritten in Fourier space as 

A similar argument to that employed in the previous subsection can be used to show that in this case Af diverges as 
at — > 0. The integral in Eq. (jBlSP gets two identical contributions from the intervals pt < and pt > 0. We will also 
take into account that since the integrand is positive, the integral is also positive and greater than the same integral 
restricted to a smaller domain of integration, so that we have 



^ I r dpte-^W^ ^ (B19) 

3a3 (19207r)(27r)4 ata^ ^'''^> 



The last integral is divergent if one takes at ^ (at least for ar ^ ^^). Thus, Af diverges unless at ^ 0. 

Using the previous result, it is easy to discuss whether a smeared version of the actual noise kernel, given by 
Eq. (|B1[) . also diverges when at 0. Each derivative in Eq. (jBl[) gives rise to an additional factor involving the 
momentum associated with the corresponding component. Additional factors involving powers of p^, p^, Py and p^ 
(odd powers of pt, Px, Py or p^ give a vanishing contribution) leave the argument employed in the previous paragraph 
unchanged and the same conclusions obtained when Cf — + hold for the smeared noise kernel as well. For example 
one has A/tut ~ ^/crtaj.. This result is in agreement with that obtained in Ref. (s^ 

On the other hand, when at the smeared noise kernel is finite even for cr^ = 0. This can be seen by taking 
(Tr = in Eq. ((BTS)) . We then have 

which is finite for at =/= 0. It is also clear that the same conclusion applies to the smeared version of the noise kernel, 
with TVtttt - 1/af. 



APPENDIX C: GENERALIZATION TO CURVED SPACE, ARBITRARY HADAMARD GAUSSIAN 

STATES AND GENERAL SMEARING FUNCTIONS 

The results obtained for the Minkowski vacuum in fiat space can be generalized to curved space and arbitrary 
Gaussian Hadamard states. They also apply to more general smearing functions. This will be shown in this appendix. 

The key ingredient is the fact that the Wightman function for any Hadamard state has the following form in a 
sufficiently small normal neighborhood of an arbitrary spacetime:^'* 

G+{x,x') ^ "(^'^') +vix,x')liia+ix,x')+w{x,x'), (CI) 
a+[x, x') 

where a+{x^x') is the geodetic interval (one half of the geodesic distance) for the geodesic connecting the pair of 
points X and x' with an additional small imaginary component added to the timelike coordinates (this prescription 



The divergence of A/" as crt —> can also be proven for ar = (i.e., in the absence of smearing along the spatial directions) by taking 
(Tr = in Eq. IIBlSt and replacing with an arbitrary but fixed positive value in Eq. IIB19I I. 

For a general spacetime it may not be guaranteed that the series has a non-vanishing radius of convergence rather than being an 
asymptotic series 1621 . However, for an analytic spacetime it can be proven that the radius of convergence is non-zero for globally 
Hadamard states WA 16311 . We will restrict ourselves to analytic spacetimes in this appendix, which is anyway the case for the spacetimes 
considered in the rest of the paper. 
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will be defined more precisely below); u, v and w are smooth functions with v and w expandable as 

oo 

v{x,x') = ^Vn{x,x')a'''{x,x'), (C2) 

n=0 

CXD 

w{x,x') — ■Wnix,x')a"{x,x'). (C3) 



n=0 



where u{x,x')^ Vn{x,x') and Wn{x,x') satisfy the Hadamard recursion relations, which uniquely determine u(x,x') 
and v{x, x'). On the other hand, wqIx, x') is not uniquely determined and contains the information on the particular 
state that one is considering, the remaining Wn{x,x') are also determined once a particular choice of woixjx') has 
been made. Note that for the Minkowski vacuum in flat space v(x, x') and w{x, x') vanish and u{x, x') is simply given 
by a constant. Other Hadamard states in flat space have a non-vanishing w{x,x') while u(x,x') and v{x,x') remain 
unchanged. 

The biscalar functions u{x,x'), Vn{x,x') and Wn{x,x') can in turn be expanded in the following way: 



m=0 



where Furthermore, it will be convenient to employ Ricmann normal coordinates, which can always be 

introduced in a normal neighborhood. Given a set of normal coordinates {y^}, the geodetic interval can be simply 
expressed as a{y,y') = (l/2)(j/^ — y'^){y'^ — y"^)'ntiv- The prescription for a+ mentioned above corresponds then to 

a+{y,y') = {l/2)[-{y^-y'^-ief + {y-y'n_ 

Given a smearing function s(x), one can consider the following smearing of the product of two Wightman functions 
N{x,x') = Y{,e[G+{x,x')G+{x,x')]: 



N= / d^x y^-gix) s{x) / d'^x' y^-gix') s{x') N{x,x'). (C5) 



Next, one changes from a set of absolute coordinates for each one of the the two points at which the kernel N{x,x') 
is evaluated to a set of absolute coordinates for the first point and a set of relative coordinates for the location of the 
second point with respect to the first one. In particular we will choose Riemann normal coordinates for the relative 
location of the second point (in order to study the divergences in the coincidence limit it is sufficient, by taking small 
enough smearing sizes, to consider small enough convex neighborhoods where both normal coordinates can be defined 
and the series in Eq. (|Gip is convergent). Eq. (|G6p then becomes 



M= / d^Tv/=5(^s(x) / d%y/^§:{y)S4y)N{x,y). (C6) 



For simplicity, in our discussion below we will consider smearing functions Sx (y) which are Gaussian and independent 
of X. However, in Sec. IC3l we will explain how our results can be extended to more general smearing functions. We 
also note that we will not analyze the integrals in x since those should be finite (we are considering globally Hadamard 
states and regular smearing functions): only the integrals in y are relevant for the UV divergences associated with 
the coincidence limit, which corresponds to y — > 0. There could still be IR divergences arising from the integrals in 
x, but we will be considering smearing functions which decay sufficiently fast so that this is not the case. 

Taking all the previous considerations into account, it becomes clear that when calculating the smeared kernel Af, 
the most divergent terms from N{x, x') = Re[G+(a;, x')G'^{x, x')] will be of the form l/a\{x, x'), lncr+(a:, x')/a+{x, x') 
or 1/ u+(x^x'). When expressed in normal coordinates, the contribution to N due to a term of the first kind, which 
will be denoted by Ni and corresponds to the product of two l/cr+(x, x') terms, coincides with the expression for the 
Minkowski vacuum in flat space. Hence, one can directly apply the results obtained in Appendix A. Furthermore, one 
expects that the leading divergence to the smeared function M in the various limits of vanishing smearing sizes will 
come entirely from that kind of terms and that the other terms will only give rise to subleading divergences. If this is 
true, the main conclusions in Appendix A will also apply to general Hadamard Gaussian states in curved spacetime. 
Let us, therefore, study more carefully the contributions to M from the different kinds of divergent terms and check 
that this is indeed the case (other possible divergent terms in addition to the three kinds mentioned above correspond 
to multiplying one of those three by some positive power of cr°, and they will be discussed in Sec. IC 1 bp . 

In order to analyze the contribution from the other two kinds of divergent terms, denoted by M2 and A/s, we will 
proceed analogously to Appendix A and work in Fourier space for the relative normal coordinates. We start by 
considering a term of the form \na^{x,x') / a+{x^x'). Using the same Gaussian smearing functions as in Appendix A 
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(as explained above, we do not include any dependence on the absolute coordinate for the first point), its contribution 
to Af when smearing around a null geodesic is given by 



(2^) 



dp,dp,4'pe-P>"e-P"-''"e-P"''" i [L{p) + (Li^p))* 



(C7) 



where L{py,pu,p'^) is the Fourier transform of ln(cr-|_(y,y')/A^)/(T-|_(j/,y') and we need the second term inside the 
square brackets because we are interested in the real part of the product of two Wightman functions. An explicit 
expression for L(pt,,p„,P^) can be determined using the following two Fourier transforms: 



a+{x,x') 

\n{(T+{y,y')/\^) = -2 / d'^p e^P'-^-y'^ lim 



d^pe'P^y-y'^9{~p")6{-p^), 



(C8) 



ln(2m-'A^) + 7- 1 



sHp) 



(C9) 



where Eq. (jC9P was derived in Appendix A. 2 of Ref. J64|- We can then write L{py,pu,p'^) as a convolution of those 
Fourier transforms: 



L{p) 



4(27r)3 / d^q9{-p" + q"))S{-{p-qf) 



X lim 



1 

TT 



) ^ 

dm? 



ln(2TO2A2) + 7- 1 



After a certain amount of calculation, one obtains the following result for L{p): 

1 



L(p)-2(27r)3|(?(-pO)p/ 



p^ 



[ln(2A2)+7-l]0(V)<5'(V) , 



(CIO) 



(Cll) 



where Vf denotes Hadmard's finite part prescription (a generalization of the principal value prescription) whose 
precise definition can be found in Refs. [6J, [gl, [6^. We will use this result in the next subsections to compute in 
Fourier space the contribution Af2 to the smeared kernel J\f for different kinds of smearings and analyze under what 
conditions it is finite. 

On the other hand, from Eq. (jCSp one can immediately see that the contribution from the term l/a^{x,x') is 
simply given by 

1 



2tt 



dpvdpud pe 



(C12) 



1. Smeciring around null geodesies 

a. Contributions from M2 and A/3 

As we found in Appendix A, when considering a smearing of the noise kernel for the Minkowski vacuum around a 
null geodesic, N diverges as l/tT„cr^CT^ in the limit of small tT„ (as long as CTj. ^ 0, otherwise it diverges as l/tr^cr^). 
Whereas the contribution Mi from terms of the form l/cr^(a;,x') will exhibit the same divergent behavior, we will 
show in this subsection that M2 and M-i , the other contributions to Af, are arbitrarily smaller than JVi for sufficiently 
small fT„ or a^. 

Let us start with Afs- First, one rewrites Eq. (jC12[) as follows: 

-^3 = 7^ [ dp,dpud^pe-P"^''"e-P"''"^e-P"''U{puPv - f) 
Ztt J 

= / dp J dpue~P>"e'P>" d\p\e-\P\"^U{,/^~\p\) 
Jo Jo Jo 

poo poo 

I I 2222 2 

= / dp. / dpue-P^^^e-P^'^^e-P-P^^^ 
Jo Jo 




(C13) 
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where we introduced the new variable ^ = PvPu in the last equality. Next, using the positivity of the integrand, the 
fact that the original integral was invariant under interchange of p„ andp„, and the fact that value of the exponentials 
is always equal or less than one, one can derive the following bound for small CTu and CTi,: 



A/'s < -T^(\n(7u + \na,,) + 0(1), 



(C14) 



where the higher-order terms involve positive powers of au and (when is not small, one has an expansion only 
in terms of cr„ and the Inert, term is absent). On the other hand, for (7^ = the bound is As < (7r^/8)(l/cr„(Tt,). 
Let us now turn our attention to Af2- Substituting Eq. (jCll[) into Eq. (jC7p . we get 



dpydpud pe 



o 2 2 2 2 ^2 2 



Vf 



0{PuPv -P ) 



PuPv - P^ 



+ [ln(2A2) +7-1] S^ipuPv (C15) 



The contribution from the second term inside the curly brackets has the same form as As. Hence, we need to 
concentrate on the first term. After a lengthy calculation one can derive the following bound: 



M2<L 



1 



0(ln(T„, In cr„), 



(C16) 



where L is a constant of order 1. On the other hand, for (7^ = the bound is M2 < L/au<Jy. We can see that M2 
dominates over A/3 in the limit of small (t„ or Uy. Nevertheless, it is still A/i that provides the leading contribution in 
that limit. 



b. Remaining terms 

The contribution from the remaining terms in N{x,x') to the smeared function J\f, which can be seen from terms 
of the different kinds that we have already analyzed multiplying them by some positive powers of a and ct" can be 
analyzed as follows. If one is working with Riemann normal coordinates and the corresponding Fourier variables, 
each a factor will give rise to a momentum d'Alembertian (l/2)ri^y{d / dp^){d / dp^) acting on the Fourier transform 
of that term in N{x,x'). Similarly, each factor will give rise to a linear differential operator {d/dp^) acting on 
the Fourier transform. Integrating by parts in the Fourier space expression for M so that the d'Alembertian acts on 
the Fourier transform of the smearing functions, it will produce a factor of the form {—a^a^PuPv + "^cTrP^ ~ 2ct^). 
Proceeding analogously with the linear operator (d/dpf^), one gets the factors 2(T^p„, 2alpv and ^.a^pj (where the 
index j corresponds to one of the two orthogonal spatial components) when /i equals u, v or j respectively. 

Next, one needs to see how the results of the integrals in Sec. IC 1 al change when the integrand is multiplied by 
positive powers of p„, p^, and pj. Odd powers of pj for any j give a vanishing result since the integrals and integrands 
there are symmetric under a sign change oipj, whereas every (pj)"^ factor can be written as {—l/2)(d/da'^). Similarly, 
every factor or pi can be written as {—l/2){d/da'^) or {~l/2){d/dal) respectively.^^ On the other hand, having 
an odd power of Pu only or p^ only gives a vanishing result because the integrals and integrands in Sec. IC 1 al are 
symmetric under a simultaneous sign change of Pu and py. However, odd powers of PuPv do not vanish in general 
because the integrands are not symmetric under interchange of pu or py only. Therefore, one needs to check how the 
main results for the integrals in Sec. lC 1 al change when the integrands are multiplied by apuPv factor. In Eq. (jC13() it 
gives rise to a factor ^, which implies an additional l/cr^ factor multiplying the final result for A/3 in Eq. (jC14p . One 
can similarly find that an additional PuPv factor in Eq. (jC15[) also implies a 1 /ct^ factor multiplying the final result 
for A'a in Eq. (ICTHl) . 

We are finally in a position to discuss the effects of a and ctq factors multiplying the contributions to the smeared 
kernel J\f. Each power of p^ and pf, (including their a* and a* accompanying factors) will typically give rise to cr^ 
and cr.y factors respectively. A PuPv factor (with its accompanying cr^cr^ factor) will give rise to a a^a'^/a'^ factor. 
And each p^ factor (with its accompanying <t^ factor) will give rise to a factor. Thus, we see that the divergent 
behavior in the limit of small (t„ remains unchanged or even gets improved. (An analogous conclusion would apply 
in the limit of small ar ) It then follows that the behavior of Af in the small CTu limit is still dominated by the flat 
space vacuum contribution A/i . 



In general when two functions satisfy a certain inequality that does not imply that their derivatives will satisfy it. This will, however, 
be the case when considering the bounds derived in the previous sections in order to analyze the leading divergent behavior in the limit 
(Tu — » 0. 
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c. Smearing of the actual noise kernel 

The actual noise kernel involves a number of functions multiplying the kernel N{x,x') — Kc[G^{x,x')G^{x,x')] 
and differential operators acting on it. When using relative Riemann normal coordinates for the location of the second 
point with respect to the first one, the part of these linear operators that depends on the relative coordinates of the 
second point can be entirely expressed in terms of functions and tensor fields (such as the metric and the curvature 
tensors) as well as partial derivatives. That is not the case in general for the operators associated with the first 
spacetime point, but the dependence on the first point does not exhibit a divergent UV behavior. 

The functions multiplying N{x, x') can be expanded in terms of the relative Riemann normal coordinates, which 
involves powers of a and cr" that can be treated as explained in the previous subsubsection. As we saw, they either 
leave the divergent behavior in the limit of small au unchanged or decrease the degree of divergence. 

On the other hand, the partial derivatives {d/dy^) simply correspond to ip^ factors in Fourier space, which can 
also be dealt with as discussed in the previous subsubsection. Even powers of p„ increase the degree of divergence in 
the limit of small a^- every factor gives rise to a 1/cr^ factor. All others leave the degree of divergence unchanged. 
Moreover, since the p\ affect both the kind of terms contributing to A/i and those contributing to N2 and A/a, 
the conclusion that the leading divergent behavior when tT„ —> is given by the Minkowski vacuum result remains 
unchanged for the actual noise kernel. 



2. Smearing along timelike geodesies and on spacelike hypersurfaces 



The results for smearings of the noise kernel on spacelike hypersurfaces and along timelike geodesies obtained for 
the vacuum state in Minkowski can be generalized proceeding analogously to what was done in the previous subsection 
for smearings along null geodesies. 

Let us start by considering As : 



AA3 = J dptd^pe-P>"e-p"''Uip^t-f) 

= 4 / dpte-p"''" / d\p!\ \p\e-\^"''" 5{pt - \p\) 
Jo Jo 



Hence, we have 



4 / dptPte-P'^'e-P'"-. (C17) 



^f3 = (CI8) 



which is finite provided that Ct 7^ or (7^ 7^ 0. 

Let us now turn our attention to J\f2. Applying the result in Eq. (jClip to this case, we get 

^f2 = -J dptd^pe-P'y^e~P'^'''\vf 0{pi - f) + [ln(2A2) + 7 - l] S\p', ~ f)^ 



(C19) 



The contribution from the second term inside the curly brackets has the same form as A3. Hence, we only need to 
concentrate on the first term. After a slightly lengthy calculation, one gets the following bound for A/'2 (when ct^ 7^ 0): 

AA, < ^ + (C20) 

where Ci and C2 are positive dimensionless constants which are finite provided that at ^ (they behave like — Incrj 
in the limit of small at)- On the other hand, for ar — one has a bound given by 

m < 4, (C21) 

where D is some positive dimensionless constant which is finite provided that at (it behaves like — In at in the limit 
of small at). This shows that just a temporal smearing is enough to render Af2 finite. Of course A/'2 will also be finite 
if, in addition to the temporal smearing, some but not all of the (orthogonal) spatial directions have a non-vanishing 
smearing size. 
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Proceeding similarly to Sees. IC 1 bllUTcl one can argue that factors involving positive powers of a and cr° as well as 
the derivative operators in configuration space do not alter the main conclusion in this subsection. Thus, a smearing 
along the timelike direction is enough to have a finite result for the smeared versions of iV(x, x') and the actual noise 
kernel. 



3. More general smearing functions 



In Sees. [CT] and [C2] we considered Gaussian smearing functions for the relative normal coordinates. However, as 
pointed out at the beginning of this appendix, when transforming from a pair of sets of absolute coordinates for the 
two points where the noise kernel is evaluated to one set of absolute coordinates and one set of relative ones, the 
form of the smearing functions will change in general. That will also happen when transforming to Riemann normal 
coordinates if one had initially chosen a different kind of relative coordinates. Furthermore, even in the flat space 
case one may be interested in considering other kinds of smearing functions. For instance, one may wish to consider 
a smearing function adapted to a spherical surface rather than a plane. 

The results obtained in these appendices can be extended to more general smearing functions, making it possible 
to cover the situations described in the previous paragraph. The essential idea is simple: the previous results can be 
generalized for smearing functions which can be locally approximated by the Gaussian smearing functions in Riemann 
normal coordinates considered so far. More specifically, by "locally approximated" we mean that in those coordinates 
the new smearing functions can be expressed as the Gaussian smearing functions times a factor involving an expansion 
in positive powers of a and tr". The procedure described in Sec. IC 1 bl can then be employed to show that the main 
results remain unchanged. Moreover, the detailed form of the smearing functions for large values of the Riemann 
normal coordinates (or even outside the range where they can be defined) is not relevant when concerned with the 
divergent behavior (when certain smearing sizes tend to zero) as we are are here. Or put in a different way, even 
though in general the form of the smearing function can be significantly distorted when transforming from the original 
coordinates to relative Riemann normal coordinates, this effect becomes less and less important when considering 
sufficiently small smearing sizes (which is in any case the relevant regime to study the UV divergent behavior): for 
sufficiently small scales the coordinate transformation is characterized by the linear map between tangent spaces (the 
Jacobian matrix evaluated in the coincidence limit) plus small corrections involving positive powers of a and ct", 
which can be dealt with following the procedure described in Sec. IC 1 bl 

These points can be illustrated with a simple example. Consider a spatial smearing function adapted to a sphere 
in flat space with the following form (in spherical coordinates): 



fir,e) 



1 



(2^)3/V,r2a2 b{ae) 



exp 



{r-rof 




" 02 - 


2^2 


exp 


. 2a2_ 



(C22) 



where b{ag) is some dimensionless function which ensures the proper normalization of the smearing function and 
tends to 1 when erg — > [note also that for 9 we actually have half a Gaussian since its domain is (0, tt)]. This can be 
written in terms of Riemann normal coordinates (cartesian coordinates in this case) adapted to the plane tangent to 
the sphere at the point (r = tq, 6* = 0) as follows: 



1 



(2^)3/2a,r2a2 b{ae) 



exp 



exp 



y 



1 + 



{z-rof {x^ +y'^){z~ro) {x^ + y^){z - ro) 



(C23) 



We can see that for any given ro > 0, if one chooses sufficiently small <Tr and ae, the higher-order terms become 
negligible (in the region where there is not a large suppression due to the exponential factors they are very small). 
Thus, f{x,y,z) corresponds to a Gaussian smearing function in cartesian coordinates with Uj, = ctj, = roag and 
= cTr plus small corrections. 
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